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Abstract 
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extension of the Standard Model (MSSM) and their implications for searches 
at accelerators. 
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I. INTRODUCTION 


In the seminal papers of Refs. it was pointed out that the Earth’s motion around 
the Sun can produce a sizeable annual modulation of the signal in experiments of direct 
search for heavy relic particles. 

Actually, the analysis of a new set of data, recently collected by the DAMA/Nal Col¬ 
laboration (in the period denoted by the Collaboration as running period # 2) supports 
the possible presence of an annual modulation effect in the counting rate for WIMPs: the 
hypothesis of presence of modulation against the hypothesis of absence of modulation is sta¬ 
tistically favoured at 98.5% C.L. Remarkable features of this measurement, obtained with 
an exposure of 14,962 kg x day, are: 

i) An analysis of the experimental data, based on a maximum likelihood method, pins 

down, at a 2-cr C.L., a well delimited region in the plane ■C'^icaiar^ where is the 
WIMP mass, crgcX*r°”^ is the WIMP-nucleon scalar elastic cross section and ^ = Px/Pi is 
the fractional amount of local WIMP density Px with respect to the total local dark matter 
density pi. This modulation region is shown in Fig. 1, which is reproduced 

here from Fig. 6 of Ref. 0 (the values of plotted in Fig. 1 are normalized to the 

value Pi = 0.3 GeV cm“^). The ensuing l-a ranges for the two quantities are: = 59^14 

GeV and x nb §. 

ii) The new data conhrm a previous indication of an annual modulation (at the 90% 
C.L.) found by the same Collaboration, by using a smaller sample of data, collected in the 
running period 1, with an exposure of 4,549 kg x day IQ. Most remarkably the 2-a C.L. 
region from data of Ref. |Q] is entirely contained inside the 90% C.L. region derived from 
data of Ref. [0, also shown in Fig. 1 (the open solid curve denotes the 90% C.L. upper 
bound derived in Ref. p, by using pulse shape analysis). 

iii) Because of the property at point ii), when the data of the two running periods (with a 

total exposure of 19,511 kg x day) are combined together, one obtains a more delimited 2-a 
C.L. region in the plane “ '^x’ which is fully embedded in the previous regions. 

Consequently, the determination of and remains very stable: nix = 591^4 GeV 

and 'C<^scriar°^^ = ^ 10”^ lit) (if pi is normalized to the value pi = 0.3 GeV cm“^). By 

combining the two sets of data, the hypothesis of presence of modulation increases to 99.6% 
G.L. 

It is noticeable that the two sets of data have been taken under different operating 
conditions, since the experimental set-up was dismounted and reassembled between the two 
running periods. 

In extracting the contour lines of Fig. 1 from the experimental data, the values of some 
astrophysical parameters (the root mean square velocity Urms of the WIMP Maxwellian 
velocity distribution in the halo, the WIMP escape velocity Vesc in the halo, the velocity Vq 
of the Sun around the galactic centre), relevant for the event rates at the detector, had to 
be chosen. The values adopted in Fig. 1 refer to the median values of these parameters in 
their experimentally allowed ranges (reported, for instance, in 0), namely: Urms = 270 Km 
s“^, Vesc = 650 Km s“^, n© = 232 Km s“^. 

In Refs. 0,^ we derived the theoretical implications of the experimental data of [Q, 
assuming that the indication of the possible annual modulation reported there was due to 
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relic neutralinos. We selected the relevant supersymmetric configurations and discussed how 
these may be investigated by indirect searches for relic WlMPs and at accelerators. 

In the present paper we apply a similar analysis to the new, much more significant set of 
data of Ref. 0 and we show that these data are fully compatible with an interpretation in 
terms of a relic neutralino as the major component of dark matter in the Universe. We pin 
down the regions of the supersymmetric parameter space relevant for this neutralino and 
derive the implications for search at accelerators. 

A word of caution is in order here. As also remarked in Ref. 0, although the new DAMA 
data appear to bring more evidence for a possible annual modulation effect, first singled out 
in Ref. this effect awaits further confirmation by additional experimental investigation 
in WIMP direct detection |^. Actually, the DAMA/Nal Collaboration has already collected 
new data over the past year; moreover, the experiment still keeps running under good 
stability conditions and is expected to provide increasingly significant statistics in the 
future. Furthermore, it is remarkable that, as subsequently discussed in the present paper, 
the supersymmetric configurations singled out by the annual modulation effect are also 
explorable at accelerators and in terms of indirect signals of relic neutralinos (i.e., in terms 
of antiprotons in space and of up-going muons at neutrino telescopes). 


II. SUPERSYMMETRIC MODEL 


In this paper we consider the neutralino as a WIMP candidate, able to induce annual 
modulation effects in direct particle dark matter searches. This supersymmetric particle is 
defined as the lowest-mass linear superposition of photino ( 7 ), zino (Z) and the two higgsino 
states (H^, H 2 ) CT 

X = ni7 T 0 , 2 ^ T (1) 


We also define a parameter P = a\ + a\ in terms of which we classify neutralinos as: 
gaugino-like when P > 0.9, mixed when 0.1 < P < 0.9 and higgsino-like when P < 0.1. 

As a theoretical framework we use the Minimal Supersymmetric extension of the Stan¬ 
dard Model (MSSM) |T^, which conveniently describes the supersymmetric phenomenology 
at the electroweak scale, without too strong theoretical assumptions. This model has been 
extensively used by a number of authors for evaluations of the neutralino relic abundance 
and detection rates (a list of references may be found, for instance, in []^). 

The MSSM is based on the same gauge group as the Standard Model, contains the 
supersymmetric extension of its particle content and two Higgs doublets Hi and H 2 . As a 
consequence, the MSSM contains three neutral Higgs fields: two of them [h, H) are scalar 
and one (A) is pseudoscalar. At the tree level the Higgs sector is specified by two independent 
parameters: the mass of one of the physical Higgs fields, which we choose to be the mass 
niA of the neutral pseudoscalar boson, and the ratio of the two vacuum expectation values, 
defined as tan/3 = {H 2 )/{Hi). Once radiative corrections are introduced, the Higgs sector 
depends also on the squark masses through loop diagrams. The radiative corrections to 
the neutral and charged Higgs bosons, employed in the present paper, are taken from Refs. 


11,12 
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The other parameters of the model are dehned in the superpotential, which contains all 
the Yukawa interactions and the Higgs-mixing term and in the soft-breaking La- 

grangian, which contains the trilinear and bilinear breaking parameters and the soft gaugino 
and scalar mass terms. 

The MSSM contains a large number of free parameters. To cast it into a form adequate 
for phenomenology, it is necessary to introduce a number of restrictive assumptions at the 
electroweak scale. The usual conditions, which are also employed here, are the following: i) 
all trilinear parameters are set to zero except those of the third family, which are unihed to a 
common value A] ii) all squarks and sleptons soft-mass parameters are taken as degenerate: 
mi, = mq^ = mo, hi) the gaugino masses are assumed to unify at Mqut, and this implies 
that the U(l) and SU(2) gaugino masses are related at the electroweak scale by Mi = 
(5/3)tan2 OwM^. 

After these conditions are applied, the supersymmetric parameter space consists of six 
independent parameters. We choose them to be: M 2 ,/i, tan/3, m^, mg, A and vary these 
parameters in the following ranges: 10 GeV < M 2 < 500 GeV, 10 GeV < |/i| < 
500 GeV, 75 GeV < m^ < 1 TeV, 100 GeV < mo < 1 TeV, —3 < A < +3, 1 < tan/3 < 
50. We remark that the values taken here as upper limits of the ranges for the dimensional 
parameters, M 2 , /i, mg, m^i, are inspired by the upper bounds which may be derived for these 
quantities in SUGRA theories, when one requires that the electroweak symmetry breaking, 
radiatively induced by the soft supersymmetry breaking, does not occur with excessive hne 
tuning (see Ref. |T^ and references quoted therein). 


Our supersymmetric parameter space is further constrained by all the experimental limits 
obtained from accelerators on supersymmetric and Higgs searches. Thus, the latest data 


from LEP2 on Higgs, neutralino, chargino and sfermion masses are used [pAjpA 

Moreover, the constraints due to the 6 —> s + 7 process (see, for instance. Refs. ig-0) 
have to be taken into account. In our analysis, the inclusive decay rate BR(i3 —> X^y) is cal¬ 
culated with corrections up to the leading order. Next-to-leading order corrections 
are included only when they can be applied in a consistent way, i.e. both to standard-model 
and to susy diagrams. This criterion limits the use of next-to-leading order corrections 
to peculiar regions of the supersymmetric parameter space, where the assumptions, under 
which the next-to-leading order susy corrections have been obtained, apply We re¬ 
quire that our theoretical evaluation for BR(i3 —X^y) is within the range: 1.96 xl0“^ < 
BR(i3 —X^y) < 4.32 xl0“^. This range is obtained by combining the experimental data 
of Refs. PB|J27[| at 95% G.L. and by adding a theoretical uncertainty of 25%, whenever the 
still incomplete next-to-leading order susy corrections cannot be applied. 

Since we are exploring here the neutralino as a stable dark matter candidate, we have 
to further constrain the parameter space by requiring that the neutralino is the Lightest 
Supersymmetric Particle (LSP), i.e. we have to exclude regions where the gluino or squarks 
or sleptons are lighter than the neutralino. We also have to disregard those regions of 
the parameter space where the neutralino relic abundance exceeds the cosmological bound, 
derivable from measurements of the age of the Universe and of the Hubble constant 

0. Gonservatively, for this cosmological bound we take < 0.7 {h is the usual Hubble 

parameter, dehned in terms of the present-day value Hq of the Hubble constant as h = 
Hq /(100 km s“^ Mpc“^)). The neutralino relic abundance is calculated here as illustrated in 
Ref. PO. Inclusion of coannihilation effects [M in the calculation of are not necessary 
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here, since the instances under which these effects might be sizeable are marginal in our 
supersymmetric parameter space. 


III. SELECTION OF SUPERSYMMETRIC CONFIGURATIONS BY THE 

ANNUAL MODULATION EFFECT 

We discuss now which region in the susy parameter space is selected by the new DAMA 
modulation data 0. 

Let us start by converting the region delimited by the 2-cr C.L. dashed contour line in 


rriy of Fig. 1 into an enlarged one, which accounts for the uncertainty 


,1 , (nucleon) 

the plane 

in the value of p/. If a possible flattening of the dark matter halo |^| and a possibly sizeable 
baryonic contribution to the galactic dark matter 
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are taken into account, the following 
range for pi has conservatively to be taken: 0.1 GeV cm“^ < Pz < 0.7 GeV cm“^. One then 
obtains from the 2-a G.L. region of Fig. 1, where the total dark matter density is normalized 
to the value pi = 0.3 GeV cm“^, the relevant 2 -ct G.L. region of Fig. 2 (hereafter denoted 
as region R). 

Now we have to find which supersymmetric configurations, out of those in the parameter 
space outlined in Sect. II, are selected by the requirement that (^x’^‘^scaiar^°^^) ^ To this 
purpose we evaluate and ^ in the MSSM scheme previously defined. 

The neutralino mass is evaluated as usual by taking the lowest mass eigenstate of the 
neutralino mass matrix 
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The neutralino-nucleon scalar cross-section is calculated with the formula 


cr, 


(nucleon) 

scalar 


j-ed 

TT 


FhR FhIr FIz Ari— I AT T-t / a 2 r-> 


mt 




( 2 ) 


where the two first terms inside the brackets refer to the diagrams with h- and iL-exchanges 
in the t-channel (the A-exchange diagram is strongly kinematically suppressed and then 
omitted here) and the third term refers to the graphs with squark-exchanges in the s- 
and u-channels 


35 . The mass m^ed is the neutralino-nucleon reduced mass and 


Fh = (—Oi sin 9w + ^2 cos 9w) (% sin a + cos a) 
Fh = (—ai sin 9w + <^2 cos 9w) (03 cos 0 — 04 sin a) 

h,H = E k^q^mq{N\qq\N). 

g 


(3) 


The angle a rotates FT® and into h and FT, and the coefficients given by 

k^-type =cosa/sin/3 and F^type = —sina/ sin/3 for the up-type quarks, and by F^.type = 
—sina/cos/3 and Fj^type = ~ cosa/cos/3 for the down-type quarks. The matrix elements 
< N\qq\N > are meant over the nucleonic state. By using the heavy quark expansion 
one may rewrite the quantity Ih,H as follows 


r(0) 
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^u—type9u “T ^d—type9d'! 


(4) 


where 
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( 5 ) 


4 19 2 23 25 

Qu = ■^{’rriN + -^o'TrAf — cict-^n), 9d = -^{rriN + + -^acr^Ar). 

Here ct^at is the so-called pion-nucleon sigma term, a-^jq = + rrid) < N\uu + dd\N >, 

and the parameter a is related to the strange-quark content of the nncleon y by 


a = y 


m.c 


niu + md 


y = 2 


< Ar|gg|Ar > 

< N\uu + dd\N > 


( 6 ) 


For these parameters we nse the following valnes: ct^at = 45 GeV y = 0.33 ± 0.09 


and 2ms/{mu + i^d) = 29 [^]; thns, nsing the central valne of y, we obtain Qu = 123 GeV 
and Qd = 288 GeV. 

In the sqnark-exchange terms of Eq. (|^) J2i denotes a snm over the mass eigenstates, 
stands for the sqnark propagators 


Pm - 2 


+ 


m 


m 


{m^ — mqY — (m^ + mq) 


(7) 


and the Aq. and Bq. coefficients are given by 


Aq^ — COs6q{Xq + Zq) + SVnOqiYq + Zq) 

Bq^ = COS 6q{Xq - Zq) + Sm9q{Zq - Yq) 

Xq = - (^cos 9wTzqa2 + sin ^ q 

type*^4 y type*^3 

type — 


K = 


sinc^H^-^ai 


7 _ 

type 


2 sin PMz 


2 cos (3Mz ’ 


( 8 ) 


where T^q, Yq^, Yq^ refer to the isospin and to the hypercharge qnantnm nnmbers of qL,m 
respectively. The conplings Aq^ and Bq^ may be obtained with the snbstitntion sin6'q —>• 


cos 9q and cos 9q 


— sin 9„ 


In onr nnmerical applications the sqnark propagators in Eq.(|3) have been regnlarized 
by inserting appropriate widths in the denominators. In general, it tnrns ont that the 
Higgs-exchange amplitndes are largely dominant over the sqnark-exchange ones, the latter 
competing with the former ones almost exclnsively when an enhancement in their size is 
originated by a mass hne-tnning in the sqnark-propagator denominators. 

As for the valnes to be assigned to the qnantity ^ = p^/pi we adopt the standard rescaling 
recipe [^]. For each point of the parameter space, we take into acconnt the relevant valne of 
the cosmological nentralino relic density. When is larger than a minimal valne 
compatible with observational data and with large-scale strnctnre calcnlations, we simply 
pnt = 1. When tnrns ont to be less than (Gh^)min, and then the nentralino may only 
provide a fractional contribntion to we take ^ The 

valne to be assigned to is somewhat arbitrary, in the range 0.01 < < 0.3. 

We nse here the valne (Gh^)inin = 0.01, which is conservatively derived from the estimate 
Ggalactic ~ 0.03. 

Using the previous formulae we hnd that a large portion of the modulation region R is 
indeed covered by supersymmetric conhgurations, compatible with all present physical con¬ 
straints. This set of susy states, which will hereafter be denoted as set S, is displayed in Fig. 
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2 with different symbols, depending on the neutralino composition. In Fig. 2(a) we notice 
that a quite sizeable portion of region R is populated by supersymmetric configurations with 
neutralino relic abundance inside the cosmologically interesting range 0.01 < < 0.7. 

Thus we obtain the first main result of our analysis, i.e. the annual modulation region, sin¬ 
gled out by the DAMA/Nal experiment, is largely compatible with a relic neutralino as the 
major component of dark matter. This is certainly the most remarkable possibility. How¬ 
ever, we also keep under consideration neutralino conhgurations with a small contribution 
to Q^h‘^ (see Fig. 2(b)), since also the detection of relic particles with these features would 
provide in itself a very noticeable information. 

The neutralino relic abundance is plotted versus the quantity terms of 

the neutralino composition in Fig. 3. Here we remark some anticorrelation between the two 
plotted quantities. This feature is expected on general grounds, as discussed for instance 
in Ref. [^. In fact, it is due to the combination of two properties: (i) the direct detection 
rate is proportional to oc , where a »nn is the neutralino-neutralino 

annihilation cross-section, (ii) usually and i functions of the supersymmetric 

model parameters, are either both increasing or both decreasing. Therefore, neutralinos 
with lower values for the relic abundance have higher couplings with matter (this feature is 
attenuated, when rescaling in is operative; this occurs here for Q^h"^ < 0.01). 

In view of the discussed anticorrelation between and it is remarkable that 

the relatively large neutralino-matter cross-sections, implied by the DAMA modulation 
effect, agree with a relic neutralino making up a major contribution to dark matter, i.e. 
with a neutralino whose relic abundance falls into the cosmologically interesting range 0.01 < 
< 0.7. Most of the neutralino conhgurations falling in this range of turn out to 
be gaugino-like. 

We further notice that recent observations and analyses point to values of hlmatter 


somewhat smaller than those considered in the past: 0.1 < hlmatter < 0.4. If we combine this 
range with the one for h\ 0.55 < h < 0.80 and require that a cold dark matter candidate 
(such as the neutralino) supplies ~ (80-90)% of hlmatter, we obtain: 0.02 < < 0.2. 

This turns out to be the most appealing interval for relic neutralinos. It is remarkable that 
this range for Q^h"^ is densely populated by conhgurations of set S (see Fig. 3). 


IV. FURTHER PROPERTIES OF THE CONFIGURATIONS SINGLED OUT BY 

THE ANNUAL MODULATION EFFECT 

Let us proceed now to an analysis of other main properties of the conhgurations of set 
S, related to a possible investigation of these supersymmetric states at accelerators. 

As is already clear from Fig. 2, the set S contains neutralino compositions of various 
nature, from higgsino-like to gaugino-like ones. This property is further displayed in Fig. 
4, where we show the location of the conhgurations of set S in the plane P.-M 2 , for two 
representative values of tan (3. 

The properties of our set S relevant to searches of neutral Higgses at accelerators are 
displayed in Fig. 5. Section (a) of this hgure shows a scatter plot of set S in term of mh 
and tan /3, section (b) provides essentially the same information, but in terms of and the 
quantity sin^(Q; —/3), which is the relevant coupling for the channels of possible neutral Higgs 
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production at LEP. In the plot of section (a) it is apparent a correlation between tan (3 and 
m/j. This is due to the fact that the rather large values of the neutralino-nucleon scalar cross- 
section, crgcaiar°”^ ~ (10“® — 10“®) nb, as required by the annual modulation data, impose 
that either the couplings are large (then large tan/3) and/or the process goes through the 
exchange of a light particle. Thus, Higgs-exchange dominance and ^ (10“® —10“®) 

nb require a very light h at small tan/3, and even put a lower bound on tan (3: tan/3 > 2.5. 


At larger values of tan /3, the mass mh is less constrained, also because, at large tan /3, 
the squark-exchange diagrams may occasionally compete with the Higgs-exchange ones in 
keeping ^ sizeable value. From Fig. 5(a) we notice that a good deal of susy 

conhgurations are explorable at LEP2, while others will require experimental investigation 
at a high luminosity Fermilab Tevatron, which should be capable to explore Higgs masses 
up to roh ~ 130 GeV P5| , ^ . 

In Fig. 6 the conhgurations of set S are shown in the plane — tan/3 denotes the 
lightest top-squark). This scatter plot reveals an interesting correlation: at small tan /3 only 
light H’s are allowed. In the Appendix it is shown that this feature occurs as a joint effect 
due to the b —>■ s + j constraint and to the annual modulation data . 


From the previous results, it then turns out that annual modulation data and 6 —s + 7 
constraint complement each other in providing stringent bounds on both ruh and at 
small tan (3. For instance, for tan 13 < 5 one has < 105 GeV and mi^ < 350 GeV. 

Finally, in Fig. 7 we display the scatter plot of set S in the plane - tan/3. Since the 
reach of LEP2 extends only up to the dashed vertical line, at ~ 50 GeV, the exploration 
of the whole interesting region will require Tevatron upgrades or LHG. Under favorable 
hypothesis, TeV33 could provide exploration up to the vertical solid line. 

Apart from exploration at accelerators, conhgurations of set S may be investigated by 
means of indirect measurements of relic neutralinos, such as cosmic-ray antiprotons il and 
neutrino huxes from Earth and Sun ( and references quoted therein). On the basis of a 
preliminary analysis, we found that conhgurations of set S provide quite signihcant signals 
in both instances. In the case of antiprotons, a large fraction of conhgurations of set S 
provide p huxes at the level of the measurements by the balloon-borne BESS experiment 
m- These conhgurations will be further investigated with the data collected during the 
Shuttle hight by the AMS experiment [^. A similar situation occurs for the neutrino huxes 
induced by conhgurations of set S, which turn out to be within the reach of MAGRO and 
Baksan 0 neutrino telescopes. Details of our analysis on the indirect detection searches 
are presented in Ref. . 

We end this section by some more general theoretical considerations. We have discussed 
here the physical implications of the annual modulation data in the framework of a MSSM at 
the electroweak scale, since this scheme provides the simplest and least-constrained model 
for discussing susy phenomenology. However, we have also performed an analysis of the 
modulation data in the framework of Supergravity (SUGRA) theories. The results of this 
study are presented in Ref. [0. We simply report here that we have ascertained that a 


fraction of conhgurations of set S are indeed compatible with SUGRA schemes, even more 
so when the unihcation conditions, which are usually imposed at GUT scale, are somewhat 
relaxed, for instance by allowing deviations from a strict unihcation assumption in the Higgs 
masses at the GUT scale [1^ . It is remarkable that these conhgurations fall into the region 






















of susy parameter space where electroweak symmetry breaking occurs without excessive hue 
tuning between competing terms. A simple case of this feature occurs for the neutralino 
mass, whose range for the annual modulation conhgurations is well within the no-£ne-tuning 
upper bound < 0(100 GeV) 0. 


V. CONCLUSIONS 

The new data of the DAMA/Nal experiment [Q, which support a possible annual modu¬ 
lation effect in the counting rates for relic WIMPs, previously reported by the same Collab¬ 
oration have been analysed here in terms of relic neutralinos. We have proved that the 
annual modulation data are largely compatible with a relic neutralino making up the major 
part of dark matter in the Universe. 

We have also investigated the possibility of exploring the supersymmetric states, selected 
by the annual modulation data, at accelerators. We have demonstrated that an analysis of 
the main features of these susy configurations is within the reach of present or planned 
experimental set-ups. In particular, we have found the following results: 

a) The sizeable neutralino-nucleon elastic cross-sections, implied by the annual modu¬ 
lation data, entail a rather stringent upper bound for in terms of tan (3. In particular, 
this property implies that no susy conhguration would be allowed for tan/3 < 2.5. A large 
portion of the region covered by the scatter plot in the plane mh - tan/3 is explorable at 
LEP2, the remaining one will be at TeV33. 

b) The annual modulation data and the 6 —»• s -|- 7 constraint complement each other in 
providing a correlation between tan /3 and the mass of the lightest top-squark. 

As remarked in the introduction, a solid conhrmation of the annual modulation effect 
as singled out by the DAMA/Nal Collaboration will require further accumulation of an 
increasingly signihcant statistics with very stable set-ups over a few years. However, it is 
worth noticing that the detection of this effect, if confirmed by further experimental evidence, 
would turn out to be a major breakthrough in establishing the existence of particle dark 
matter in the Universe. It is very rewarding that the features of this dark matter particle are 
widely compatible with those expected for the neutralino, both in MSSM and in SUGRA 
schemes, and that several of its properties can be explored in the near future at accelerators 
and by indirect searches for relic neutralinos. 
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VII. APPENDIX 


Here we discuss the origin of the correlation between rrti^ and tan (3 which is apparent in 
the plot of Fig. 6 at small tan/?. 

Let us start by considering how the 6 ^ s + 7 constraint ||I 6 |-^ correlates the three 
parameters tan/3,m/i and mi^. Thus, leaving momentarily aside the annual modulation 
data, let us consider in the plane m^^-tan/? the regions of our parameter space which satisfy 
all accelerator constraints (including 6 —s + 7 ) and the further requirement that rrih is 
below some arbitrarily hxed value mj/. In Fig. A.l these regions are represented by the 
domains on the left of the various lines, which are denoted by the following values of mj/: 

= 80, 90,100,110 GeV. It is possible to show that the b —>■ s +7 constraint is instrumental 
in establishing the peculiar shape of the various contour lines at hxed mh- 

If we now combine the plot of Fig. A.l with the one of Fig. 5(a) we obtain the situation 
displayed in Fig. A.2, the allowed region being the one on the left of the various curves, 
depending on the values of From this hgure we see how the tan/d-mj-^ correlation, 
occurring in Fig. 6 , is due to the joint effect of 6 —»• s + 7 and annual modulation data. 
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FIGURE CAPTIONS 


Figure 1 - Annual modulation regions singled out by the DAMA/Nal experiments in 
the plane The dotted contour line denotes the 90% C.L. region deduced from 

the data of the running period A 1 [0) fh® solid contour line delimits the 2-cr C.L. region 
deduced from the data of the running period A 2 [Q, and the dashed contour line delimits 
the 2 -ct C.L. region, obtained by combining together the data of the two running periods. 
The solid open curve denotes the 90% C.L. upper bound, obtained in Ref. [^, where a pulse 
shape analysis of the events was used. This hgure is reproduced from Fig. 6 of Ref. 0]. 

Figure 2 - Scatter plot of set S in the plane The dashed contour line 

delimits the 2-a C.L. region, obtained by the DAMA/Nal Collaboration, by combining 
together the data of the two running periods of the annual modulation experiment [Q]. The 
solid contour line is obtained from the dashed line, which refers to the value pi = 0.3 
GeV cm“^, by accounting for the uncertainty range of pi, as explained in Sect. Ill (the 
region delimited by the solid line is denoted as region R in the text). Displayed in this 
hgure are only the representative points of the susy parameter space, dehned in Sect. II, 
which fall inside the region R. Dots, crosses and circles denote neutralino compositions 
according to the classihcation given in Sect. 11. Sections (a) and (b) refer to conhgurations 
with 0.01 < — OT and with < 0.01, respectively. 

Figure 3 - Scatter plot of set S in the plane Dots, crosses and circles 

denote neutralino compositions according to the classihcation given in Sect. 11. The two 
vertical solid lines delimit the D^h^-range of cosmological interest. The two dashed lines 
delimit the most appealing interval for as suggested by the most recent observational 

data. 

Figure 4 - Scatter plot of set S in the plane p-M 2 . Sections (a) and (b) refer to two 
representative values of tan/3: tan/3 = 8 and tan/3 = 30, respectively. The solid curves 
denote the iso-mass curves which delimit the annual modulation region R, i.e. the iso-mass 
curves for = 34 GeV and = 107 GeV. The dashed curves denote the neutralino 
composition, and correspond to P = 0.1, 0.5, 0.9. The hatched region is excluded by LEP 
at x/i = 183 GeV. 

Figure 5 - Section (a) - Scatter plot for set S in the plane rrih - tan/3. The hatched 
region on the right is excluded by theory. The hatched region on the left is excluded by 
present LEP data at a/s = 183 GeV. The dotted and the dashed curves denote the reach of 
LEP2 at energies ^/s = 192 GeV and ^/s = 200 GeV, respectively. The solid line represents 
the 95% G.L. bound reachable at LEP2, in case of non discovery of a neutral Higgs boson. 
Section (b) - Scatter plot for set S in the plane mh - sin^(a — /3). The hatched region on 
the left is excluded by LEP data at yU = 183 GeV. 

Figure 6 - Scatter plot for set S in the plane mp - tan I3. The hatched region is excluded 
by LEP data (without any restriction on other masses). 

Figure 7 - Scatter plot for set S in the plane - tan /3. The hatched region on the 
left is excluded by present LEP data. The dashed and the solid vertical lines denote the 
reach of LEP2 and TeV33, respectively. 


11 


Figure A.l - Regions of the parameter space defined in Sect. II, which satisfy all 
accelerator constraints (including 6 —s + 7 ) and the further requirement that rUh is below 
some arbitrarily hxed value m\. The various lines denote the following representative values 
of m*y^\ m’^ = 80, 90,100,110 GeV. The allowed regions are given by the domains on the left 
of the various curves for each value of 

Figure A.2 - Allowed region in the plane - tan/? when the plot of Fig. A.l is 
combined with the one of Fig. 5(a). The hatched region on the left is excluded by LEP data 
(without any restriction on other masses). 
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